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Abstract 

According to statistical turbulence theory, the ensemble averaged squared vorticity 
Pe is expected to grow not faster than dpE/dt ~ p^ . Solving a variational problem 
for maximal bulk enstrophy (E) growth, velocity fields were found for which the 
growth rate is as large as dE/dt ~ E 3 . Using numerical simulations with well re- 
solved small scales and a quasi-Lagrangian advection to track fluid subvolumes with 
rapidly growing vorticity, we study spatially resolved statistics of vorticity growth. 
We find that the volume ensemble averaged growth bound is satisfied locally to a re- 
markable degree of accuracy. Elements with dE/dt ~ E s can also be identified, but 
their growth tends to be replaced by the ensemble-averaged law when the intensities 
become too large. 
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1 Introduction 



Fluid turbulence may be characterized as a tangle of intermittent vortices em- 
bedded in regions of straining motion [1,2]. Theoretical and computational 
studies of the evolution of isolated intense vortices [3,4], pairs of vortices 
[5,6,7,8], highly-symmetric vortex tangles or ensembles of randomly distributed 
vortices [9,10] have contributed considerably to our understanding of the sta- 
tistical properties of homogeneous isotropic Navier-Stokes turbulence, but nu- 
merical computation of high-vorticity events in a turbulent flow remains elu- 
sive because of the high demands on spatial and temporal resolution. The 
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properties of such events, their frequency and maximal intensities, are im- 
portant for small scale mixing, the efficiency of combustion processes, and for 
modeling turbulence. The situation for the inviscid Euler equation is not much 
different. Recently, Bustamente and Kerr [11] discussed in detail the sensitivity 
of vorticity growth on grid resolution and de-aliasing techniques for interact- 
ing highly-symmetric anti-parallel vortices. They came to different conclusions 
from Hou and co-workers [12] who observed a depletion of the vortex stretch- 
ing in their Fourier smoothing method. The question of finite-time blow-up in 
solutions of the Euler equations remains an open area of investigation. 

For viscous fluids there is a bound on the relation between the growth rate 
of the volume integrated squared vorticity E, namely dE/dt < aE 3 with 
known prefactor a, and there is also a calculation of optimal fields for which 
dE/dt ~ E 3 . On the other hand, the ensemble averaged squared vorticity pE is 
not expected to grow faster than dp^/di ~ p 3 ^ 2 . In an effort to understand the 
relation between these two results and how they translate to turbulence, we 
conducted high-resolution direct numerical simulations with a special focus on 
the dynamics on small scales and the evolution of strong vorticity amplification 
elements. This is the objective of the present work. 

The strong temporal variations near "almost singular" events in turbulence 
can only be resolved with a sufficiently small time step, which together with a 
high spatial resolution requires the storage of huge data files. We use a cubic 
box with periodic boundaries in all directions and solve the Navier-Stokes 
equations numerically, 



where / is a large-scaling forcing. We apply the pseudospectral method with 
2/3 de-aliasing and obtain a homogeneous isotropic and statistically stationary 
three-dimensional flow at a Taylor microscale Reynolds number of R\ = 107 
[13]. The grid size is 2048 3 , so that when expressed in terms of the Kolmogorov 
length i]k = v 3 ^ / (e) 1 / 4 (with the mean energy dissipation rate (e)) the sides of 
the box are 68377^- long and there are 3 grid points per tjk- Since the crossover 
from viscous to inertial range occurs near a scale of 8t]k, we can resolve singular 



events into the viscous range. The time step At = 0.003r^, where r v = J vj (e) 
is the Kolmogorov time, is well within the limits discussed by Donzis and 
Sreenivasan [14]. The events we would like to study are followed for a time 
interval of about 55 t v units or 1830 output steps. This would add up to about 
5 x 10 13 velocity field values that have to be stored. In order to avoid this, 
we turn to the so-called quasi-Lagrangian method [15] which eliminates the 
large scale sweeping motion superimposed on the localized singular events we 
want to study. Specifically, we follow 100 Cartesian cubes Vl simultaneously 



du 

~dt 



+ {u ■ V)u = - Vp + uAu + / 
V • u = 0, 



(1) 
(2) 




2 



through different regions of the evolving Navier-Stokes flow. The motion of 
the subvolumes Vi is fixed by the advection of a Lagrangian tracer in their 
center, and their sides are kept aligned with the outer coordinates. The boxes 
have a side length of L = 17t]k, corresponding to 51 3 grid points. This reduces 
the number of velocity field values that have to be stored to 1.4 x 1CT 3 of the 
original estimate. 



2 Analytic predictions on the growth rate of enstrophy 



The vorticity is the curl of the velocity field, uj = V x u, and the enstrophy 
is the volume integral of its intensity, 

E(t) = J u 2 dV. (3) 
v 



It follows from the Navier-Stokes equation that in any incompressible viscous 
Newtonain fluid the growth rate of the enstrophy, dE/dt, obeys ([16,17]), 



dE(t) 
dt 



ui ■ Vu) ■ udV - 2v / (VufdV . (4) 



v 



From this it can be shown that E(t) cannot grow faster than 
dt lw 6 



with v the kinematic viscosity of the fluid and c = y2/n (for details see 
[16,17]). This holds for the volume integrated quantity and does not make any 
assumptions on the flow. Incompressible flow fields that maximize the enstro- 
phy production, and thus the growth rate of enstrophy, were recently found by 
solving an optimization problem [16,17]. At high Reynolds number the max- 
imum enstrophy growth rate (5) is realized by colliding, axially symmetric 
vortex rings. At lower Reynolds numbers, maximum enstrophy generation is 
realized by interacting Burgers vortices with dE/dt ~ E 1 ^. 

There is a second analytical result that pertains to the growth rate of the 
ensemble averaged squared vorticity (or enstrophy density) 

Pe = (u 2 ) (6) 
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for the particular case of homogeneous and isotropic (box) turbulence. A direct 
consequence of the von Karman-Howarth (KH) equation [18] for the velocity 
correlations, when the volume average (•) in (7) agrees with the ensemble 
average that appears in the KH equation, is derived in [19,20] and states that 



where S is the skewness of the longitudinal velocity derivative, and u x is the 
x-component of the turbulent velocity field. It is an empirical fact that S < 0. 
It has been observed that the skewness S is basically constant for Taylor 
microscale Reynolds numbers R\ < 200 and it grows weakly as \S\ ~ -Ra' 11 
for R\ > 200 [21,22]. Thus we will assume for purposes of discussion and 
data analysis that dps/dt ~ p% 2 holds approximately. This exponent is much 
smaller than the one in the upper bound (5). 

In order to make the relation between the two relations more explicit, we 
rewrite the first one assuming that E = L 3 pE = L 3 (u) 2 ) for a box of length L 
and (e) = vpE to bring in the Kolmogorov length. Then the bound suggests 
that 

d 27c 3 ( L\ & 3/2 u . 



The key difference then is a ratio of lengths: if the sidelength L of the vol- 
ume V is of the order of the Kolmogorov scale, both results imply the same 
exponent, despite the different boundary conditions and derivations. But if 
the volume is larger than r]K, a stronger variation is possible: in a situation 
where the vorticity content of the volume is below the mean that enters the 
definition of the Kolmogorov length, the local dissipation length is larger than 
the statistical average, and a built up of entstrophy will reduce the local value, 
thereby increasing the contribution from the factor (L/tjk)- 



3 Results 



3.1 Local quasi- Lagrangian analysis 



Configurations as highly symmetric as the colliding vortex rings that realize 
the maximum instantaneous enstrophy generation cannot be generically ex- 
pected in a turbulent flow. In particular, the high-amplitude vorticity events in 
turbulence arise preferentially in the form of tubes [1,2,23] which are rapidly 
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Fig. 1. (Color online) Vortex tubes and dissipation sheets in homogeneous isotropic 
turbulence. Isosurfaces of the vorticity magnitude square (or local enstrophy) 
U = uj 2 (cyan) and the energy dissipation rate e = 2vSijSij (red) with 
Sij = (dui/dxj + duj/dxi)/2 the rate of strain tensor. Both surfaces are shown 
at the level of ten times their mean. The displayed volume is 1/16 of the full simu- 
lation box. 

stretched and deformed by background and self-induced straining motions 
[24,25]. Strong strain (or shear) results in by high-amplitude energy dissi- 
pation rates. The spatial distribution of high-amplitude events of the local 
enstrophy and energy dissipation rate, illustrated in Figure 1, underlines this 
behavior. Isosurfaces of the local enstrophy, Q(x,t) = u 2 (x,t), (cyan) at ten 
times the mean value show primarily elongated structures. The isosurfaces of 
energy dissipation rate, e(x,t) = {u /2)(dui/dxj + duj/dxi) 2 , (red) also at ten 
times the mean value, reveal sheet-like structures between the high-vorticity 
events. This illustrates that while the ensemble averaged values of energy dis- 
sipation and vorticity are related by (e) = up E this does not apply to their 
instantaneous and local values. It does show, however, that extreme events 
occur at neighboring locations. 

To capture this quantitively, we study the probability distribution functions 
(pdf 's) of local energy dissipation and vorticity in Figure 2. The pdf 's of e(x, t) 
and Q(x,t) show stretched exponential tails indicative of strong small-scale 
intermittency (see Figure 2(a)). The tail is more extended for Q than for e, in 
agreement with the observations in [26]. Fat tails imply that large amplitude 
events are significantly more probable than for a Gaussian distributed signal 
with the same second moment. In Figure 2(a) we also mark the global max- 
ima in e and Q identified within the advected volumes Vj, for three particular 
Lagrangian trajectories. The locations far out in the tails document that our 
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Fig. 2. (Color online) Statistics of local enstrophy and energy dissipation rate, (a) 
Probability density functions of local enstrophy, O, and energy dissipation rate, e, 
given in units of their means, respectively. The vertical arrows mark the global 
maxima of f2 and e, respectively, in case of three Lagrangian tracers, no. 14, 
33 and 45. (b) Joint probability density function of local enstrophy and energy 
dissipation rate. The distribution is normalized by both single quantity distribu- 
tions, p(zi, Z2) /[p{z{)p{z2)\, in order to highlight the statistical correlations between 
Z\ = e/(e) and Z2 = Q/pe- Color coding is in decadic logarithm. 



quasi-Lagrangian tracking is able to detect high- amplitude events and that 
extreme events in both quantities are spatially correlated and located within 
our advected volume. The local correlation between high-amplitude local en- 
strophy and energy dissipation events is further supported by the joint pdf 
in Fig. 2(b), where p(e/(e), Q/ pe)/[p(€/ (e))p(Q/ pe)] is shown. The maximum 
values appear in the upper right of the support where the largest amplitudes 
for both are present. High- amplitude fluctuations in energy dissipation and 
local enstrophy density are thus strongly statistically correlated and found 
very close together, in both, space and time. 



We now turn to the study of the time evolution of extreme events within our 
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Fig. 3. (Color online) Quasi-Lagrangian analysis of enstrophy. (a) Time series of 
Ei{t) for all 100 subvolumes Vl are plotted. The time traces that reach the largest 
local maxima for El are colored differently and their labels are indicated in the 
legend. Enstrophy is given in units of u 2 / L with L = ldrjK- (b) Enstrophy growth 
rate versus enstrophy. The enstrophy growth rate, d£x/di, is given in units of 
u 3 /L 3 . The a priori upper bound dEL/dt = 27\/2/(8v 3 Vir 3 ~)E L is indicated as a 
red line. The growth that follows from the von Karman-Howarth equation [18], 
dEL/dt ph — 7S/3y/lEE^ 2 with a derivative skewness of S = —0.5, is indicated as 
a solid blue line. The dashed blue line has the same slope and serves as a guide to 
the eye. Color coding is as in panel (a). 

subvolumes. Figure 3(a) shows time traces of the local enstrophy 

E L (t) = J u 2 dV. (9) 

v L 



Since we are interested in the relation between large values of d£x/d£ with 
El, we show in Fig. 3b the same data as a scatter plot in the plane spanned 
by El and dEL/dt on a double logarithmic scale. The collection of the indi- 
vidual growth histories in the subvolumes (which can vary strongly from one 
to another) is bounded from above by the scaling dEL/dt < E^J 2 , indicated 
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there by the dashed line. This shows that the local growth rate and enstrophy 
are related very much the volume averages, Eq. (4), so that the effects 

of sweeping are averaged out. It suggest that the estimate 



dE 



L 



dt 

v L 



J (w • Vu) • ujdV ~ {u 2 } 3/2 V L ~ y/El/V L . (10) 



also holds locally. We have also verified this scaling in a volume Vl/8, i.e. in a 
cube with half the sidelength. The striking observation is that the envelope of 
the local enstrophy growth follows the scaling of (7). However, we wish to stress 
that (generically) there are nonvanishing enstrophy fluxes, f 9v uj 2 u ■ dA ^ 0, 
across the bounds of Vl- Therefore (10) is only a heuristic estimate. 



3.2 Local analysis of extreme events 



The previous analysis shows that when the average volume has a diameter of 
about the Kolmogorov length or above, the bulk estimate for the enstrophy 
growth is recovered. Stronger growths, therefore, may only occur on smaller 
scales. To probe for this, we turn to the study of the time evolution of local 
extrema within the cell. We select the grid point x* with the fastest local 
growth rate within the subvolume, dQ/dt\ max = max xG v L (t)[d£l/di\. As an 
indication of the numerical uncertainty, we also show the 26 growth rates for 
points on a 3 3 cube surrounding x*. Figures 4(a) and (b) demonstrate that 
the maximum position always gives the outer envelope of the curve. The curve 
is continuous by construction, but its slope is discontinuous when the position 
x* of the point of maximal growth rate jumps discontinuously within the cell. 
Panel (a) shows also that long periods of low variability and small growth 
rates are interrupted by short violent outbursts of the local enstrophy growth 
rate. The two subsequent intervals I and II for tracer No. 45 mark exactly 
such a rapid growth event. 

Figure 4(b) shows the same data as (a) in the plane spanned by the local en- 
strophy and its growth rate, dfl/dt and Q. The first part follows a scaling that 
resembles the global enstrophy growth bound, dQ/dt ~ Q 3 . It is connected 
with a rapid stretching by two approaching and interacting tubular vortex 
segments and thereby is similar to the vortex ring collision seen for the opti- 
mal growth events. However, when dissipation kicks in and the high-amplitude 
dissipation region between the tubes forms, the growth is weaker as can be 
seen Fig. 5 where we plot a sequence of isosurface plots of Q(x, t) and e(x, t). 1 

1 Although we tracked 100 subvolumes simultaneuously within the evolving turbu- 
lent flow, this specific event was the only one observed that displayed such a rapid 
growth in connection with a large local enstrophy amplitude. 
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Fig. 4. (Color online) Local analysis of very rapid enstrophy growth events, (a) Time 
evolution of the maximum of dO/dt in Vl for tracer No. 45. Growth rates at the 
grid point of the maximum (blue for I and black for II) and the 27 neighbouring 
points (cyan for I and red for II) are shown. The curve is piecewise continuous since 
one local maximum in Vl takes over a former at a different grid point, (b) Replot of 
the data from (a) in the dTt/dt-Vl plane. The arrows in panel (b) indicate the time 
evolution. 



We might interpret it as a precursor of an "attempted" finite-time blowup. 
Further unbounded growth is prohibited by the imperfect collision of the tube 
segments, since it is highly unlikely that fully developed turbulence generates 
highly symmetric flow configurations that pin the local vorticity maximum 
fixed in space [9]. 



4 Conclusions. 



Our results show that the local enstrophy growth dE^/dt ~ E L expected 
from the volume and ensemble averaged Navier-Stokes equations can also be 
observed locally in volumes of a diameter equal to 17 Kolmogorov lengths that 
are advected by the flow. The typical enstrophy growth time scales in these 
local volumes are of the order of the Kolmogorov time scale t v . They show 
that there is negligible influence from the sweeping motion, so that enstrophy 
growth is dominated by local events. Within the boxes larger exponents can 
be observed when the total vorticity inside the volume is weak, since then the 
instantaneous, local vorticity gives a larger (local) Kolmogorov length, and the 
increase in enstrophy and corresponding decrease in the (local) Kolmogorov 
length contributes to an increasing ratio (L/t]k). During such events the lo- 
cal vorticity can be amplified with a rate dE/dt ~ E 3 , but this growth is 
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Fig. 5. Structures in an extreme enstrophy growth event. Upper six panels: Iso- 
surfaces of (red: 20 (CI), gray: 10 (fi)). Lower six panels: Isosurfaces plots of e 
(yellow: 2 (e), cyan: 20 (e)). Panels (A) for t/r v = 37.44, (B) for t/r v = 37.8, (C) 
for t/Tn = 38.16 correspond with time interval I in Fig. 4; (D) for £/t„ = 38.52, (E) 
for t/r v = 38.88 and and (F) for t/r v = 39.24 with time interval II. 

cut off after some t v . A further growth of enstrophy is then limited by the 
shear layer and thus the enhanced dissipation that is formed between the in- 
teracting vortex filaments. This conclusion is supported both by displaying 
the evolution of vortex filaments (as seen in Fig. 5) and by monitoring the 
energy dissipation rate and enstrophy accumulated in the comoving volumes 
(not shown). Moreover, they seem to be rare and overwhelmed by the E 3 / 2 
behavior. Nevertheless, already for the present moderate Reynolds number we 
found such rapid growth events and expect their more frequent appearance 
for larger ones. 
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